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Abstract 

We discuss the validity of replacing the complicated three-body confinement operator of the 



> 

. Y string junction type by three kinds of approximation which are numerically much simpler to 

o\ ; 

\ handle: a one-body operator with the junction point at the centre of mass, a two-body operator 

<^ ' corresponding to half the perimeter of the triangle formed by the three particles, and the average 



of both. Two different approaches for testing the quality of the approximations are proposed: a 



Qh! geometrical treatment based on the comparison of the potential energy strengths for the various 

1^ ' inter quark distances, and a dynamical treatment based on the comparison of the corresponding 



effective string tensions using a hyperspherical approach. Both procedures give very similar results. 
It is shown how to simulate the genuine string junction operator by the approximations proposed 
above. Exact three-body calculations are presented in order to compare quantitatively the various 
approximations and to confirm our analysis. 
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I. THE Y STRING POTENTIAL AND POSSIBLE APPROXIMATIONS 



The QCD lattice calculations support the idea that the confining potential energy in a 
three quark system can be, at best, simulated by the so-called Y-shape potential In this 
scheme, the three pointlike quarks (located at the apices of a triangle ABC), are connected 
by three flux tubes starting from a junction point (J) in such a manner that it minimizes the 
sum of the distances to the three quark positions. Note that effective QCD theories support 
also this idea Q|. 

If the angle corresponding to an apex is greater than 120°, the junction point is precisely 
at this apex, whereas if all the angles are less than 120°, then the equilibrium junction 
position coincides with the so-called Toricelli (or Fermat or Steiner) point; it corresponds to 
a point for which the corresponding angles AIB, BIC, and AIC are all equal to 120°. 

Thus the genuine string potential in a baryon is deflned by 

Vy = cr{AB + AC) if A> 120° {+ circular permutation on B and C), 

Vy = a{IA + IB + IC) if A,B,d< 120°, (1) 

where a is the string tension, in principle a constant of the theory. This expression has been 
used in several works (see for instance Refs. 0,0,0,0)- The success of simple quark models 
using a potential to simulate the conflnement without gluonic modes can be justified by the 
fact the lowest gluonic excitation energy in the three quark system is found to be about 
1 GeV in the hadronic scale jlo| . 

Efficient methods to deal with Y-shape interaction relies either on Monte-Carlo algorithms 
(see for instance Refs. [ijIsD or hyperspherical methods (see for instance Ref. js]). However 
it is very difficult to implement in a numerical code, essentially for two reasons: i) it is a 
three-body operator which makes the equations very complicated; ii) the presence of two 
different expressions depending upon the configuration makes the integration domain not 
evident to handle. This complexity explains why, for practical calculations, the genuine Y- 
shape potential is often replaced by other expressions that are considered as approximants. 

The most popular one is to consider half the perimeter of the triangle formed by the three 
quark coordinates and to define the confining potential as 

Vc = ^iAB + BC + CA), (2) 
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which is also called the A-shape potential. However it is shown in Ref. that the three 
quark potential energy is better represented by the Y-shape interaction than by the A-shape 
one. Nevertheless, potential is a two-body operator, free from the complication due to 
angles: hence it is much more simple and it is still widely used in practice (see for instance 



Refs. 



A Aj ■ Xj colour dependence associated with a two-body linear confinement does produce 
such a 1/2 factor. Although this colour prescription is perfectly relevant in the case of one- 
gluon exchange, there is no theoretical justification to apply it for the confinement potential. 
This 1/2 factor is close to the value 0.53, predicted by lattice calculation if one tries to 
replace the Y-shape by a A-shape. 

It is very important to stress that both the genuine string potential (0) and the two-body 
confining potential Q are of geometrical essence, depending only on the position of the 
quarks and independent of their masses. 

Another approximation, that was suggested, is to replace the true junction point / by the 
centre of mass G of the three quark system ^ . In this case the corresponding confining 
operator is simply 

VG = a{GA + GB + GC). (3) 

This approximation is particularly interesting since this potential is a one-body operator 
free from angle complications; hence its numerical treatment is quite easy. In contrast to 
the previous expressions, this approximation does depend on the system via the centre of 
mass coordinate. 

In this paper, assuming that Vy represents the true physics, we want to study the quality 
of the approximations (j2I), ©and their relevance. To this end, we propose two approaches: 
a geometrical one and a dynamical one relying on the hypercentral formalism. With this last 
technique, it is possible to obtain directly an average value of the confining potential energy 
depending only on one lenght parameter. Then we check their validity with an exact three- 
body treatment based on a complete hyperspherical treatment (beyond the hypercentral 
approximation). A simulation of three-body potential by sums of two-body potentials was 
performed in Ref. Q|, but with a philosophy very different from the one developed here. 

In the next section, the geometrical approach is presented. In Sec. IIIH the hyperspherical 
formalism is used to calculate the value of the effective string tension for a given system. 
Sec. IIVI is devoted to a three-body treatment of the confining potential and the simulation 



of the genuine string operator. Conclusions are drawn in the last section. 



II. GEOMETRICAL APPROACH 
A. Configuration of the system 

We are interested in the ratio of the potential energy for two forms of confining interactions 
as function of the configuration of the system. Since this ratio is scale independent and 
since the dynamical constant a disappears in this ratio, it is always possible to rescale the 
quark triangle putting AB — 1 and to deal only with the remaining apex of the triangle. 
Denoting by L the minimal distance from the junction point /, by D the distance from the 
centre of mass G, and by P half the perimeter of the triangle, we will study the ratios: 



In order to obtain analytical expressions (already complicated!) for these ratios we restrict 
ourselves to a system with 2 identical quarks of mass m located in A and B and a third one 
of mass M — xm located in C. The region < x < 1 corresponds to QQq systems with a 
light and two heavy quarks; the region x > 1 corresponds to qqQ systems with a heavy and 
two light quarks. The case x = 1 corresponds to systems qqq with three identical masses. 

Since AB is fixed, the only freedom for the geometry of the system is the position of the 
apex C, that can be defined by two angles a = A and b = B. Following the previous remarks, 
one must calculate L(a, 6), P(a, 6), D{a, h, x). In fact this study is still too complicated. One 
can simplify it a lot noting that there is a symmetry versus the mediating line of AB, with the 
consequence that i?(a, h) — R{b, a). This implies that an extremum of those ratios, the only 
important thing for our consideration, always lies in the mediating line and it is sufficient 
to restrict the study to isosceles triangles. Thus, we only compute Ry/c{o,), Ra/c{o,,x) and 
RG/Y{a,x). 

In a quantum mechanical treatment of the confinement, the wave function of course 
explores all the configurations for the triangle so that, may be, the most important physical 
quantity to be computed is the average of the previous ratios. One defines 



Ry/c = Vy/Vc = L/P, Rg/c = Vg/Vc = D/P, Rg/y = Vg/Vy = D/L. 




(4) 
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B. Various distances 



The calculation of the various distances results from geometrical properties in a triangle 
and presents no difficulty. It is easy to show that 



The genuine string distance is given by 

1 



L(a) 



cos a 



if < a < 7r/6, 



L{a) = ^(tan(a) + V3) if 7r/6 < a < 7r/2. 



(5a) 

(5b) 



Half the perimeter is equal to 



P(a) 



1 + cos(a) 
2 cos(a) 



Centre of mass string distance is given by 

tan(a) + A/x2tan^(a) + (2 + x)^ 



D(a, x) 



2 + x 



(6) 



(7) 



C. genuine potential to half perimeter 

With the expressions © and (jHl), one has -Ry/c(a) = L{a)/ P(a). This ratio is always 
larger than 1, meaning that the genuine string junction potential is always more repulsive 
than the sum of the two-body confining potentials. However the ratio is 1 for a fiat triangle 
or for an infinitely stretched triangle, and presents a maximuiiij 2/^3 ^ 1.155, for an 
equilateral triangle. This is in agreement with the result of Ref. jpy|]. 

The average value (jH) in this case is equal to 



2 2 



1 - + 2 In 



1.086. 



(8) 



The value 
from QCD 



^2 Ry/c ~ 0.54 must be compared with the corresponding value 0.53 derived 
jj . Thus the average error replacing the string tension operator by the sum of the 
two-body confining potentials is of the order of 8%; this approximation can be considered 
as a good one. 
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D. Centre of mass junction approximation to half perimeter 



There does not exist special angle conditions in this case and, following formulas ^ and 
(jni), the ratio is given by RG/c{ci,x) = D{a,x)/P{a). 

The maximum value is obtained for an infinitely stretched triangle and for an infinite 
mass asymmetry. In this special case, the ratio equals 2, but in the physical part of the 
domain, this ratio is much closer to unity. A more reliable estimation results from averaging 
following the procedure The integral is cumbersome but can be evaluated analytically 
using for example the Mathematica package: 

Rc/cix) = ^(^2\x) { 2VTTx arctan (2^/l + x)/x^ 

-(2 + x)E(X)+a;2K(X)/(2 + a;)}, (9a) 

where 

X - f^. (9b) 

(2 + x)^ 

and where KfX) and E(X) are the complete elliptic integrals respectively of first and of 
second kinds jl7|. The function Rg/c{x) is presented in Fig. ^ It is always greater than 
1, indicating that the centre of mass string always overestimates the sum of the two-body 
confining potentials. For x = {QQq systems) the value is 1.168, for x = 1 {qqq systems) it 
is 1.136, and for x = 20 {qqQ systems) it is 1.248. Curiously it passes through a minimum 
1.132 for X ~ 0.600 corresponding to ssu and ssd systems. Replacing the two-body potential 
by the centre of mass string induces an error of about 15-20%. 



E. Centre of mass junction approximation to genuine potential 

In a similar way, we define the ratio RG/Yi^jX) = D{a,x)/L{a). A maximum value 2 is 
obtained in very extreme situation a ^ n /2, x —>■ oo,but, in general, the values of expression 
this ratio are very close to 1. 

The integration of expressions Ry/c{o) over the angle a is very cumbersome, the final 
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result looks like 

Rg/y{x) = ;^^^|^(^ + 3v^ln(2)) + + ^ E(vr/6,X) 

- ^Vl + x + x^ \n{Y) + Vl + x ( ^ ln(T) + ^ arctan(f/) ] j> (10a) 



')} 



where X is given by Eq. ()9b|) and the new quantities Y, T, and U are equal to 

y _ x(2Vl + x + x'^ - y/3x) 



V^{1 + x) + ^/{l + x + x^){3 + 3x 



X 



2^ 



y ^ V3 + 3x + x2- v/3(l + x) 
V3 + 3x + x2 + ^3(1 +a;)' 

t/ = ■ (lOd) 

The behaviour of this ratio is shown on Fig. |21 Some remarkable values are 1.075, 1.048, 
1.149 for X = 0, 1, 20. The function Rg/y{x) passes through a minimum 1.043 for x ~ 0.585, 
corresponding again to the S baryon. It is very instructive to emphasize that for a large 
domain of masse ratios < a; < 5, the error introduced by replacing the Toricelli point by 
the centre of mass is less than 10%. 



F. Another approximation 

Replacing the genuine string junction potential by a sum of the two-body confining poten- 
tials is a rather good approximation in any case. Replacing it by a sum of one-body centre 
of mass string potentials is even a better approximation for equal quark mass systems and 
one light-two heavy quark ones. This approximation becomes slightly worse (although not 
dramatically) for one heavy-two light quark systems. In most cases, these approximations 
are better than 10%. 

It is interesting to remark that the genuine string tension is always comprised between 
half perimeter and centre of mass junction. This last property is obvious since the Toricelli 
point is precisely the one which minimizes the sum of the distances; in contrast the former 
property is by no means obvious. One can take benefit of this remark and to define a new 
ratio as 

Rm/y{x) = ^ {Rc/Y + Rg/y{x)) . (11) 
8 



In Eq. (fTTj). the function Rg/y{x) has been computed before. The value of Rc/y can be 
computed analytically. The result is 



Rc/Y = -[^ + {l + V?>) 



TT 



2 ln(2) + ln( - 1) - 3 ln( + 1)1) ^ 0.923. (12) 



The curve Rm/y{x) is plotted on Fig. El One can remark that the values of Rm/y{x) differ 
from unity by less than 3% for all relevant values of the x parameter. So this procedure 
to simulate the genuine string junction potential seems preferable to the previous discussed 



ones. 



III. HYPERSPHERICAL APPROACH 
A. hyperspherical coordinates 

The hyperspherical formalism is an economical way to tackle the three-body problem. 
We refer to speciahzed papers for technical aspects (see for instance Ref. 0|). Here we just 
recall what is needed for our purpose. 

Let us define a reference mass m and introduce the dimensionless quantities Ui = mi/m, 
Uij = uji + ujj and u = Ui + uj2 + ujs- The first thing to do is to introduce the Jacobi 
coordinates 

pij = aij{ri -Tj), Xij = pij ( ^ ^^^^ -rk) , (13) 



(jjij 



with 



,UJiUJj UJkUJij UJ1UJ2UJ3 , . 

^ii = \ ^' Pij = \ ^' ^ = aijPij = \ • (14) 



The normalization quantities have been set in order to obtain nice properties under particle 
permutations. From now on, we particularize to the 1-2 pair and drop the 12 index every- 
where, so that the Jacobi coordinates for our problem are noted simply p (instead of P12) 
and A (instead of A12). 

Each inter-distance r^j = — Vj can be expressed only in terms of p and A, so that 
half the perimeter of the quark triangle is expressed as P{p^, A^, p ■ A). The same is true 
for the positions relative to the centre of mass Si = Vi — Rem-, so that the centre of mass 
string distance is D{p^, X^, p ■ A). Finally, the same property is valid for the genuine string 
distance L{p^, p ■ A). Over the six original variables defining the configuration, three 
have disappeared corresponding to the three Euler angles giving the orientation of the plane 



of the quarks in a fixed reference frame. The confining potential is expressed in terms of 
two distances p and A, and one angle x = (p, A) between p and A. Instead of p and A, 
the hyper spherical formalism introduces the hyperradius R and the hyperangle 6 through a 
polar transformation 

p = Rsm9, X = Rcos9 with Q<9<tt/2. (15) 
The hyperradius is invariant under quark permutations 

R = ^^^T^ =^ph + AI3 = VpIi + Ail. (16) 
The elementary volume element with the hyperspherical coordinates is simply 

dv = R^dRdn^^\ dn^^^ = cos^esin^edo sm xdxdn^^\ (i?) 

where dil^^^ is the volume element on the hyperangles and dQ^^^ the usual volume element 
on Euler angles. One has obviously 

An important property of the hyperspherical formalism is that a dominant part of the 
interaction comes from the hypercentral approximation of the potential V{p, A) which is the 
average of the potential over the hyperangles 

V{R) = ^ J V{p,\)dn^'l (19) 

For potentials invariant under rotations, as the confining term we arc interesting in, the 
expression does not depend on Euler angles so that, in practice, the hypercentral potential 
is simply 

Q /"'r/2 PIT 

V{R) = - / cos^ sin^ 0de V{R, 9, x) sin x dx- (20) 
TT Jo Jo 

In the hypercentral approximation the complicated three-body problem reduces to a differ- 
ential equation including V{R). Within this approximation, we obtain directly an average 
value of the potential confining energy depending only on the hyperradius. Such a treatment 
is feasible for the Y-shape as well as its approximants. 

In this paper we focus only on the genuine confining potential Vy and its two approxima- 
tions Vc and Vg- In term of the hyperradius it! and hyperangles 9 and %, they are defined 
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by 



VY{R^e,x) = aL{R,e,x), 

Vc{R,9,x) = aP{R,9,x), 
Vg{R,9,x) = aD{R,e,x)- 



(21a) 
(21b) 
(21c) 



The hypercentral approximations of formulas ()2H) resulting from equation ()20j). are denoted 
here VV(i?), Vc{R) and Vg{R). 

The expressions for P and D are easy to obtain. For the moment, let us quote the 
following formulas 

" ^ (22) 



D{R,9,X) = Y.—^^^ 

i<3 



and 



(23) 



The expression for L is much more cumbersome. As was mentioned above, there are three 
special cases for which the angles 9j = (f^, f^) are greater than 120° (— 1 < cos 6j < —1/2); 
in this case L = rij + rj^. The "normal case" (all Bj less than 120°) is L = rji + r/2 + tj^ 
and is less easily obtained. Let us quote the final expression: 



First region 



-1 < 



— tan + cos Y 



^2 



tan^ e + 2— tan 6* cos y + 1 



<-2' 



Li{R,e,x) = Rcose 



+ — 4 / ^ tan^ 9 + 2 — tan fc* cos x + 1 



a 



Second region 



-1 < 



— tan 6 — cos y 

UJ2 



tan^ 6-2— tan 61 cos y + 1 



<-2' 



L2{R,9,x) = RcosO 



tanO 1 /fi2 Q 

h — W — 9 tan^ — 2 — tan fj cos x + 1 



Third region 



-1 < 



1 + ( ^ - ^ ) tan ^ cos X - — 



tan2 6' 



(24) 



(25) 



(26) 



(27) 



= < --, (28) 

l + 2^tan^cosx+ gtan^^Jl -2£tan^cosx+ gtan^^ 2 
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Rcos 9 



1 + 2 — tan cos x H 5- tan 6 



1 — 2 — tan 9 cos x H j tan 6* 



(29) 



Fourth region or "normal region" 



-l/2<cosei<l Vi, 



(30) 



L4{R,0,x) = RcosO 



tan^ 6 + 



UiUJ2{ujI — UJ2) 



LU2 ~ UJl /— 

COS x + V 3 sin X 

UJ12 



tan 6' 1 , , 



Since each distance is proportional to the hyperradius, all the confining potentials under 
consideration have the form (/ = C, G, Y) 

Vi{R) = bi{uJi,uj2,U3)R. (32) 

One sees that in the hyperspherical formalism the confining potential remains linear, with 
the string tension which is no longer a but a modified value b which depends on the system. 
This dependence is not of geometrical essence, as in the previous Section for Vq, but of 
dynamical character that comes from the choice of the Jacobi coordinates. The main effort 
of this section is devoted to the calculation of this new "string" tension b. 

To stick to the previous Section and also to get analytical results, we restrict from now on 
to the case of two identical particles (particles 1 and 2 with the reference mass m = mi = 1712 
and a particle 3 with mass rris = xm). In this special case, the coefficients a and (3 in Eqs. 
(HH) are 

a = ^, p=J-^, n = aP= .r^. (33) 
x/2 \2 + x' ^ y2 + x ^ ' 

Let us study separately each approximation. 
B. Centre of mass junction 



Among the three contributions to the potential ()22|) . one is particularly simple. It corre- 
sponds to S3 = {[3Rcos6) /0J3. The averaging is trivial and gives the value ?>2[3R/ {Ibnuj^). 
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The calculation of Si is longer but presents no special difficulty. A trick to fasten the answer 
is to remark that R is an invariant under quark permutations. Consequently, the averaging 
can be performed with the hyperspherical angles 6'23 and X23 as well, so that the result follows 
immediately from the previous one giving the contribution 32/?23-R/(15vrco'i). Similarly the 
contribution due to S2 is 32/53ii?/(157rco'2). Replacing the quantities l3ij,ujk by their values 
()33|) gives the final result for the effective string tension in this approximation 



C. Half the perimeter 

The procedure is essentially the same using now potential ()23|) . The contribution due 
to = p/a is very easy to obtain; the result is 32i?/(157ra). Switching to the permuted 
hyperangles gives immediately the contributions due to r23 and r^i, namely 32i?/(157ra23) 
and 32_R/(157ra3i) respectively. Replacing the quantities aij by their values gives the 
final result for the effective string tension in this approximation 



D. Genuine string junction 

Because of the special cases, several zones of the plane (6*, x) must be isolated correspond- 
ing to the conditions ((j21I), (jSHI), (j2Hl)), where the integrand has a special form ((j2Sl), (I2ZI), 
((2ni))- These various zones are separated from the "normal zone" by a line which corresponds 
to the limit cosGj = — 1/2. It is convenient to take 6 as abscissa and x as ordinate in the 
plane. 

The condition ()24|) writes explicitly 



This curve starts with the value 27r/3 for 6 = and ends up to vr for = = arctan (1/^^). 
In the ffist region, above xi, the expression L = Li applies. 




(34) 




(35) 



Xi{0) = arccos 



Sfitan^ - V4 - 3fi2 tan^ 9 
4 



(36) 
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The condition 



X2{e) 



writes explicitly 

3(] tan ^ + yi^^Wtm?^ 



arccos 



Ti-xi{e). 



(37) 



This curve starts with the value vr/S for ^ = and ends up to for 6 = 6q. In the second 
region, below X2, the expression L = L2 P7j) applies. 

Finally a third region is found from condition ()28|1 . which writes explicitly in terms of 
two functions 



X3 (0) = arccos 



v/(3^^2tan2^- l)(3-(]2tan2 

2n tan e 



Xt{0) = 7r-xs{e). (38) 



These curves start with the values and vr for 9 = 9q and join to the common value 7r/2 for 
9 = 9i = arctan(v^/f^). In the third region, at the right of Xa? the expression L = L3 (I29|l 
applies. 

In the rest of the plane, the "normal" expression L4 is vahd. The situation is 
summarized on Fig. 

To obtain 6y(x) we integrate the function L{R, 9, x) in the plane 9, x- The result is given 
as the sum of seven contributions 



^ J LiR,9,x)dn^'^ 



8aR 



TT 



(Jl + I2 + I3 + h + h + h + h) ■ 



(39) 



We were not able to obtain an analytical formula for /. The best that we can do is to 
transform the double integral into a single integral noting that the integrals over x can be 
calculated analytically. Instead of integration over 9, it is much simpler to perform the 
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change of variable Z — Qtan^. With these options, one has exphcitly 



Z'^dZ 



rx2{z) 

I 2Z+ x/l + Z2 - 2ZC0SX sinx^x, 

rxi{Z) r I -I 

/ \/l + 3Z2 + 2V3Zsinx 

JyoiZ) . 



h 
h 
h 



(Z2 + Q2)7/2 

' Z^dZ 

{z' + n'y/'j^,^z) 

^ Z^dZ 

X3"(^) r 



sinx^^X, 



/ Uz + v/l + ^2 + 2Zcosxl sin x ^x, 



(40a) 
(40b) 
(40c) 



X 



j ^Vl + Z2 + 2Zcosx + Vl + ^^ -2ZcosxJ sinx(ix, 

sinxc^X, 



73 



Z2 rfZ 



/ \/l + 3Z2 + 2v^Zsinx 



X 



1 (Z2 + fi2)7/2 
1 (Z2 + l]2)7/2 

Vl + ^2 + 2^cosx + + -2Zcosx 
Z'^dZ 

(Z2TJ1T^ 



(40d) 
(40e) 



xtiz) 



sinxdx, 



(40f) 



X 



Vl + ^2 + 2Zcosx+ Vl + ^2 _ 2Z 



cosx 



sinxcix- 



(40g) 



Due to the symmetries of the problem, it is easy to show that Ji = Is and I 4 = Iq. The 
integrals /i, I4, If can be calculated analytically, but not I2 and I^. 

The calculation is not trivial, but some help is provided with the Mathematica package. 
One obtains 

, , , 32cT ^ f 28 + 39Q2 + 18Q^ 

Oy(x] = r^v2 • 



IStt 



4(4 + 3Q2)2 



+ 



120 + 368Q2 + 612Q4 + 641^*^ + 377^^8 + 105^^° + 9^^^ 

2Q(4 + 3Q2)2(3 + Q2)2(l + Q2)3/2 

^Vl^ dZ Z{1 + VsZ) 



AVsJo (^^ + 1^2)7/2 



X 



(l + 3Z^)E(^-f,. 



(1 - VSZfF 



X2 

2 ' 



u 



dZ Z{l + ^/3Z) 



4^/3 A {Z^ + 1^2)7/2 



X 



1 + 3Z^)E ( ^ - ^, J - (1 - V3Z)^F r^-%u 



(41a) 
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with 



X 



2 + ,T' 



uiZ] 



X2{Z) = arccos 



3Z + V4^^3Z2 
4 



{1 + Vszy' 

XaiZ) = arcsin 



2Z 



(41b) 



and F(0, m) and E(0, m) are respectively the elhptic integrals of first and second kinds 



E. Comparison of the effective string tensions 

It is interesting to compare the string tensions bcix) bc{x) (jHSI), byix) (jHj) for the 

three expressions of the confining operator. Contrary to the previous geometrical approach 
they all depend on the system; as in the preceding section we introduce the ratios of the 
string tensions ry/cix) = bY{x)/bc{x), rc/cix) = bG{x)/bc{x), rG/vix) = bG{x)/bY{x). 



1. Genuine to perimeter 

The ratio rY/c{x) is presented on Fig. |S| It is always greater than one, indicating that 
the genuine tension is more repulsive than the two-body approximation. It starts from 1 for 
X = 0, increases to a maximum 1.099 for x ~ 1 and then decreases slowly around 1.085 for 
large values of x. Those values are very close to the constant value 1.086 obtained in the 
geometrical approach. 



2. Centre of mass to perimeter 

The ratio rG/c{x) is plotted on Fig. Q It starts from 1 for x = 0, increases rapidly 
up to 1.15 for X ~ 0.5 and then presents a plateau at this value around x ~ 1; it tends 
asymptotically to 2\/2/{l + \/2) ^ 1.1715 for an infinite value of x. Although the forms of 
the curves for the hyperspherical formalism and the geometrical approach are not identical, 
the values for the corresponding ratios are rather close. 
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3. Centre of mass to genuine 



The ratio rG/vi^) is shown on Fig.|21 It starts from 1 for x = 0, raises to a local maximum 
1.0556 for X ~ 0.195, then passes through a local minimum 1.051 for x ~ 1 and raises very 
slowly to 1.08 for large values of x. Here again the values are in nice agreement with the 
geometrical approach, and are even closer to the genuine case. 



4-. Another approximation 

In the hyperspherical formalism, we find the same features concerning the approximations 
than in the geometrical approach: for a given value of the string constant a, the centre of 
mass junction potential is more repulsive than the potential with the genuine junction which 
is in turn more repulsive than half the perimeter potential. These two approximations can 
be considered as quite reasonable, differing by no more than 10% as compared to the genuine 
string junction one. 

As in the case of the geometrical approach, one remarks that 6y(x) is always comprised 
between hc{x) and hG^x). Thus it is tempting to define a new approximation huix) and a 
new ratio rM/Y{,x) similar to that of Eq. (fTT|l by 

I . bnix) {bcjx) + bG{x))/2 \ , . 

The curve rjvf/y(x) is plotted on Fig. 01 One can remark from this Figure that rM/rix) 
differ from unity by less than 2% for all values of the parameter x. Like in the geometrical 
case, using this prescription to simulate the genuine string junction potential seems more 
preferable to the ones discussed previously. 



IV. THREE-BODY CALCULATIONS 



A. Strategy 

For a system composed of two identical particles of mass m and a third particle of mass 
M = xm, the three-body equation in the hypercentral approximation to the hyperspherical 
formalism, is very simple; it looks like 

[K{x, R) + b{x)R]^{R) = E^{R). (43) 
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In this expression E is the energy eigenvalue, and K{x, R) symbohzes the differential op- 
erator corresponding to the kinetic energy plus all the hypercentral potentials except the 
confining hypercentral interaction which is explicitly written as b{x)R, as it was proved in the 
preceding section. Up to now, we considered three types of confining potentials Vj, labelled 
by an index I {I = Y for the genuine three-body confinement, I = C for the two-body half 
perimeter confinement and I = G for the one-body centre of mass junction). The important 
point is that, for these three possibilities, the operator K{x, R) is the same. For a given 
string tension a, the only difference lies in a different value of the constant b{x) appearing in 
Eq. ()4H|1 . For a genuine string tension, the value b{x) = 6y(x) given by Eq. ()4H1 . while for 
centre of mass junction and half perimeter confinement one must employ bcix) and bc{x) 
of formulas and ()H5|1 respectively. The numerical results of the various possibilities are 
of course different although they should be close, within less than 10%. This property was 
the conclusion of our two previous sections. 

At the light of the behaviour of the curves Rm/y{x) and rM/y(a;), it is natural to try to 
simulate the potential Vy by defining a new confining potential 

yM = l {Vc + Vg) , (44) 

with Vc and Vq given respectively by formulas Q and ©. Let us emphasize that the unique 
above definition agrees with both the geometrical approach for the ratio Rm/y = Vm/Vy 
and the hyperspherical approach using buix). 

As the ratios Rm/y{x) and rM/yix) are very close to 1, one can hope that the Vm and the 
Vy will give similar spectra with the same string tension. The situation is not so favourable 
for potentials Vc and Vq since the corresponding ratios R and r can differ from 1 by more 
than 10%. Nevertheless, we can try to simulate results obtained from the Y-shape potential 
by using a renormalized value of the string tension in potentials Vc and Vc- 

Indeed, let us suppose that we perform an hyperspherical calculation for / approximation 
with a string constant aj. Now let us perform an hyperspherical calculation for J approx- 
imation not with an identical string constant aj = aj, but with a string constant modified 
in the following way <7j{x) = <Jirj/j{x) (crj(x) = ai bi{x)/bj{x)). Obviously, one recovers 
the original equation and thus the original results. In other words, one can simulate the 
results of a treatment based on an approximation J, by performing a treatment based on an 
approximation J, provided we change the string constant in a consistent way. This conclu- 
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sion is perfectly exact for the hypercentral approximation. This does not mean that it must 
remain exact if we perform a more sophisticated three-body treatment. The quahty of this 
simulation, as well as the results coming from Vm, are the subjects of this section. 

Our numerical algorithm to solve the three-body problem with potentials Vc, Vg and Vm 
is based on an exp ansion of the wave function in terms of harmonic oscillator functions with 
different sizes [l8|. It was checked with other methods and was proved to give results of good 
accuracy if the expansion is pushed sufficiently far (let say up to 16-20 quanta). Moreover it 
can deal easily with a relativistic kinetic energy operator. The detail of technical aspects is 
not the subject of this paper and can be found elsewhere For the present purpose, it is 
enough to say that we are able to solve in a very fast and precise way a three-body calculation 
either with a non-relativistic or relativistic expression for the kinetic energy operator. 

One-body and two-body operators are easy to implement, but three-body operators are 
much more complicated to handle, specially if one must distinguish several integration do- 
mains, as it is the case for the Y-shape potential. In this paper, the three-body problem with 
potentials Vy is solved by the hyperspherical method without the limitation of the hyper- 
central approximation f^. At the present stage, only S-wave states can be computed with 
good accuracy. This is why the simulation of a three-body operator either with a two-body, 
a one-body operator, or a mixing of both is important. 

We consider a system composed of 3 quarks of type n (for u or d), s, c, b. In this paper, 
we are only interested in the possibility to simulate the genuine confinement by simpler 
potentials and, for that purpose, it is enough to restrict the interaction to the confining 
term only. Forgetting about Coulomb, hyperfine and constant potentials, our results cannot 
be compared to physical systems. But the comparison of the various simulations between 
themselves are very instructive. But just to have some connections with the real systems, 
we put arbitrarily the masses of the quarks at the physical values (in GeV) [11]: m„ = 
0.330, rris = 0.550, rric = 1.850, nib = 5.200. The only potential taken into account is 
the linear confining potential as defined in the first section. To see the sensitivity of the 
results versus the kinematics, we perform two types of calculation: one based on a non- 
relativistic expression (Schrodinger equation), and one based on a relativistic expression 
(spinless Salpeter equation). To test also the sensitivity to excited states, we performed the 
calculations not only for ground states (L = and = 1), but also for the first radial 
excited state ( L = and = 2), and when it is possible the first orbital excited state 
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(L = 1 and = 1). 



B. Comparison centre of mass to perimeter 

This section does not deal with the Y shape potential, but nevertheless it is important 
because we compare here two treatments that we are able to handle easily, in any physical 
situation, and with a good accuracy. The possibility to simulate one treatment by the other 
is very instructive and can be tested carefully, so that firm conclusions can be drawn. 

For a number of systems, exploring a large domain of the x parameter, we calculate first 
the baryon binding energies obtained with the centre of mass junction [G approximation) 
and with the half perimeter confinement {C approximation) using the same value of the 
string tension a = 0.2 GeV^. This value is close to the accepted value coming from lattice 
calculations Then we keep this value of cr, do the calculations for the G approximation 
using a modified value of the string tension based on the arguments of the geometrical 
approach cr/j(x) = a Rc/ci^) and on the arguments of the hyperspherical approach ar{x) = 
crrc/G{x)y and compare with the C calculation with the string tension at the value a. We 
also perform the reciprocal calculations. Our quantitative results are presented in Table H] 
for a nonrelativistic kinematics. 

Note that the definition of the ratio r//j (see sec. IIII E|l implies that rj/i = l/rj/j. This 
property is not exact for the ratio Ri/j. However, the value I/Rq/y ~ 1.083 (coming from 
formula (O) must be compared with Ry/c ~ 1-086 (formula (jH))). These values are very 
close and this gives us confidence to use Rjjj ~ l/Rj/j for all cases. 

If we compare column (1) and column (4) from Table U we can see that, with the same 
string tension, masses obtained by potentials Vc or Vq are rather different. The binding 
energy can differ by about 200 MeV or more. So it is relevant to answer the question: "Is it 
possible to simulate each potential by the other, simply in adjusting the value of the string 
tension?" . 

If we look at columns (1) and (2) for systems characterized by x = 1, we can see that 
all masses are identical (actually, they differ at the 6th digit, which is at the limit of the 
accuracy of our calculation method). It is then possible to simulate perfectly the potential 
Vc with a string tension a by the potential Vg with the string tension orcjc- This property, 
which is exact at the hypercentral approximation, is also verified for the full calculation for 
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symmetrical systems. The situation is less favourable for systems with x ^ 1. In these cases, 
we can remark that the masses of column (2) are always greater than the masses of column 
(1). For snn and nss baryons, the agreement is good, especially for the ground state. But 
for very asymmetrical systems such as bnn and nhh baryons, there can exist greater mass 
differences, up to about 100 MeV. Let us remark that, for the ground state, the agreement 
is still reasonable (better than 2%). 

If we compare now data from columns (1) and (3), we can see that, in general, the 
potential Vc is not so well simulated by the potential Vq using Rc/g instead of tc/g- In 
particular, masses from these two columns are not identical for systems with x = 1 (since 
an exact simulation is impossible with the geometrical approach, we consider that using 
1/Rg/g instead of Rq/c in column (6) do not spoil our conclusions). In most cases, masses 
of columns (3) are greater than masses of column (2). But for very asymmetrical systems, a 
peculiar mass in column (3) can be closer to the one in column (1) than the corresponding 
mass in column (2). 

We retrieve the same features in comparing data from columns (4) with data from columns 
(5) and (6). But one can remark that the masses of column (5) are always smaller of the 
masses of column (4), when they are not identical. Moreover the masses of column (6) are 
generally smaller than those of column (5). This situation is the opposite of the one for 
columns (1), (2) and (3) 

Some calculations have also been performed with a relativistic kinematics, but the same 
conclusions can be drawn. We have just remarked that the agreement between masses is 
slightly less good than in the nonrelativistic case. To be complete, let us mention that, with 
this kinematics, for systems with x — 1, the differences between masses computed with Vg 
(Vg) and those computed with Vq (Vg) and the string tension multiplied by tg/g {''^g/g) 
appear at the 5th digit, which can be considered as relevant for the accuracy of our method. 

The simulation based on the renormalized string constant originated from the hyperspher- 
ical formalism gives generally good results, better than those coming from the geometrical 
approach. Nevertheless, for very asymmetrical systems, better results can be obtained with 
this last method. 

It is important to emphasize that using renormalized string tensions provides, in any 
case, much better results (6% in the worse case) than keeping a single value of a (17% in the 
worse case, 10% in the best). This study gives a strong confidence to use also a renormahzed 
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a to simulate the genuine string junction. 
C. Simulation of the genuine junction 

We consider the genuine confining three-body potential with a constant string tension 
a and its various approximations. Our quantitative results are presented in Table |n] for a 
nonrelativistic kinematics. 

In this paper, the spectra of the potential Vy is obtained by a hyperspherical formalism 
containing grand momenta K = 0,2,4. For the moment, only binding energies of S-wave 
states have been computed. They are reported in column (1) and are used as a reference to 
test the quality of the various approximations. All other values have been computed with 
a harmonic oscillator basis up to 20 quanta. The accuracy of all binding energies is better 
than 1%. 

In column (2), binding energies of the potential Vm are presented. They are obtained 
with a value of string tension which is the same that the one used for Vy- One can see 
that the numbers of columns (1) and (2) differ generally by less than 20 MeV, with only a 
notable exception, the first excited S-wave state of the bnn system. We cannot say nothing 
about the P-wave states, but we can expect that the results from the Vm potential are also 
close to those of the Vy interaction. Thus, it appears that the potential mixing equally the 
half perimeter and the centre of mass junction simulates quite well the genuine Y-shape 
interaction. 

Let us now discuss the quality of the spectra obtained with Vc and Vg and with a 

renormalized string tension as explained previously. We performed the calculation with a 

(c) 

two-body confining potential of type C and with a string constant either ar (x) = ary/c{x) 
(column (3)) or ctj^ = a Ry/c (column (5)). Then we redo the calculation with a one- 
body confining potential of type G and with a string constant either ar (x) = o"ry/c<(x) 
(column (4)) or a^^\x) = a Ry/dx) (column (6)). 

By comparing data from columns (3) and (4) in Table m we can see that the masses for 
baryons nnn and bbb are the practically the same (differences are at the level of the 6th digit 
as described in the previous section) and very close to the value of column (1). This means 
that, for this kind of symmetrical systems, the simulation of the potential Vy by interactions 
Vc and Vg, based on hyperspherical formalism, gives very good results. The situation is 
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less favourable for asymmetrical systems. For snn and nss baryons, the agreement between 
columns (3) and (4) is still good, especially for the ground state. But for very asymmetrical 
systems such as bnn and nbb baryons, there can exist greater mass differences, up to about 
100 MeV. Again, for the ground state, the agreement is still reasonable. Let us remark that 
masses in column (4) are always greater than the corresponding ones in column (3). The 
values of these two columns are in reasonable agreement with the reference results of column 



Within the geometrical approach, the results of the two procedures of simulation are 
never in perfect agreement, as we can see in comparing columns (5) and (6), but the values 
obtained generally enclose the value of column (1). Let us note that the differences between 
a mass in column (5) and the corresponding one in column (6) are generally of the same 
order than the gap between values of columns (3) and (4). 

The renormalization of the string tension as suggested by the geometrical and the hyper- 
spherical formalisms allows to compute with potentials Vc and Vq binding energies which 
are close to the ones obtained with the potential Vy- It gives in any case much better results 
than the ones obtained keeping a fixed value of the string tension. 

To be complete, it is worth mentioning that the same conclusions can be obtained with 
a relativistic kinematics, although in this case we cannot compute the reference energies 
corresponding to the potential Vy- 

D. Particles with different masses 

Up to now, the analysis has been done for at least two identical particles. General 
analytical formulas in the case of three different particles are not available, except bo and 
be in the hyperspherical formalism. In this case, we noted that the tensions are expressed 



A first possibihty to treat easily the problem of three different masses is to use the 
hypercentral approximation with bu — {be + bo)/'^-! with the above expressions for be and 



(1)- 



by 




(46) 



(45) 
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ba- This procedure would avoid a double numerical integration to get by, whereas allowing 
results better than 2%. 

Another possibility relies on the Vm potential proposed in the previous section. We have 
verified that this approximation works well in the case of two or three identical particles. 
We can reasonably assume that it will be also good in the case of three different masses. 

V. CONCLUSION 

In this paper, we studied three approximations of the genuine three-body confinement: a 
two-body potential Vc equal to half perimeter of the triangle formed by the three particles, 
a one-body potential Vq with the junction point at the centre of mass, and a mixing of 
both Vm — {Vc + Va)/2. Two approaches were investigated to test the quality of these 
approximations: a geometrical one for which the important quantities are the various dis- 
tances in the plane of the particles, another one based on the hyperspherical formalism. 
Both give very similar and consistent conclusions. The potential energy Vq overestimates 
the potential energy of the genuine junction by about 5% in most cases, and about 10% 
in extreme asymmetrical situation. The confining potential energy Vc underestimates the 
potential energy of the genuine junction by about 8%. Keeping the same value of the string 
tension in approximants can induced a 100 MeV error in the calculated masses, as compared 
to the spectra obtained from the Y-shape confinement. In this respect, the Vm interaction 
simulates the potential energy of the Y-shape interaction to better than 2%. 

Thus, using Vm with the same string tension than the genuine junction gives very good 
results at the level of the spectra. To obtain similar quality (sometimes a bit better or a 
bit worse) for the potentials Vc and Vc, it is necessary to renormalize the string tension by 
a mass dependent factor that can be analytically computed in the cases of two and three 
identical particles. This is very important to simplify the technical effort. 
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TABLE I: Simulation of interaction Vc by interaction Vq, and vice versa, for a nonrelativistic 
kinematics. Binding energy in GeV of various baryons as a function of the total orbital angular 
momentum L, the total angular momentum and parity J^, the principal quantum number N, the 
type of confinement potential (Conf.), and the value of the string junction a in GeV^. The total 

spin is equal to J and the total isospin is the lowest one. The mass ratio x is also given for each 
system. Data columns are numbered to make the discussion easier. 



system 


L 




N Conf. 


Vc 


Va 


Vg 


Vg 


Vc 


Vc 








a 


0.2 


0.2rc/G 


0.2/Rg/c 


0.2 


0.2rG/c 


0.2 Rg/c 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


nnn 





1/2+ 


1 


1.912 


1.912 


1.933 


2.104 


2.104 


2.081 


(x = l) 






2 


2.633 


2.633 


2.662 


2.898 


2.898 


2.867 




1 


1/2- 


1 


2.332 


2.332 


2.358 


2.567 


2.567 


2.539 


bbb 





3/2+ 


1 


0.763 


0.763 


0.771 


0.839 


0.839 


0.830 


{x = l) 






2 


1.050 


1.050 


1.062 


1.156 


1.156 


1.143 




1 


1/2- 


1 


0.930 


0.930 


0.940 


1.024 


1.024 


1.013 


snn 





1/2+ 


1 


1.819 


1.821 


1.827 


2.004 


2.002 


1.996 


{x = 1.667) 






2 


2.480 


2.505 


2.513 


2.758 


2.730 


2.721 




1 


1/2- 


1 


2.192 


2.214 


2.221 


2.437 


2.413 


2.405 


hnn 





1/2+ 


1 


1.652 


1.673 


1.604 


1.854 


1.831 


1.909 


{x = 15.76) 






2 


2.171 


2.300 


2.205 


2.548 


2.405 


2.509 




1 


1/2- 


1 


1.939 


2.037 


1.953 


2.256 


2.149 


2.241 


nss 





1/2+ 


1 


1.719 


1.721 


1.744 


1.894 


1.891 


1.866 


{x = 0.6) 






2 


2.335 


2.359 


2.390 


2.595 


2.569 


2.536 




1 


1/2- 


1 


2.067 


2.086 


2.114 


2.296 


2.275 


2.245 


nbb 





1/2+ 


1 


1.247 


1.278 


1.246 


1.373 


1.339 


1.374 


{x = 0.063) 






2 


1.514 


1.599 


1.559 


1.718 


1.626 


1.668 




1 


1/2- 


1 


1.397 


1.465 


1.428 


1.574 


1.501 


1.540 
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TABLE II: Comparison between genuine potential Vy with a constant string tension a, the con- 
fining potential Vm with the same string tension, and the Vc and Vg confining interactions with 
renormalized string tension (see text). Binding energies are obtained for a nonrelativistic kinemat- 
ics. Same kind of data as in Tabled 



system 


L 




Conf. 


Vy 


Vm 


Vc 


Vg 


Vc 


Vc 








a 


0.2 


0.2 


0.2 ry/c 


0.2 ry/G 


0.2 Ry/c 


0.2/RG/y 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


nnn 





1/2+ 


1 


2.032 


2.009 


2.036 


2.036 


2.020 


2.040 


(x = l) 






2 


2.785 


2.768 


2.804 


2.804 


2.782 


2.810 




1 


1/2- 


1 




2.451 


2.483 


2.483 


2.464 


2.488 


bbb 





3/2+ 


1 


0.810 


0.801 


0.812 


0.812 


0.806 


0.814 


(x = l) 






2 


1.111 


1.104 


1.118 


1.118 


1.110 


1.121 




1 


1/2- 


1 




0.978 


0.990 


0.990 


0.983 


0.993 


snn 





1/2+ 


1 


1.933 


1.913 


1.935 


1.937 


1.921 


1.928 


{x = 1.667) 






2 


2.626 


2.625 


2.639 


2.666 


2.620 


2.653 




1 


1/2- 


1 




2.317 


2.332 


2.356 


2.316 


2.345 


bnn 





1/2+ 


1 


1.752 


1.760 


1.747 


1.769 


1.745 


1.695 


(x = 15.76) 






2 


2.310 


2.390 


2.295 


2.432 


2.293 


2.330 




1 


1/2- 


1 




2.110 


2.050 


2.153 


2.049 


2.063 


nss 





1/2+ 


1 


1.826 


1.807 


1.828 


1.831 


1.815 


1.841 


{x = 0.6) 






2 


2.474 


2.468 


2.484 


2.509 


2.467 


2.523 




1 


1/2- 


1 




2.183 


2.199 


2.219 


2.184 


2.231 


nbb 





1/2+ 


1 


1.313 


1.312 


1.296 


1.329 


1.317 


1.316 


{x = 0.063) 






2 


1.645 


1.618 


1.574 


1.662 


1.599 


1.647 




1 


1/2- 


1 




1.488 


1.453 


1.523 


1.476 


1.509 
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FIG. 1: Ratio Rc/ci^) from the geometrical treatment (see formula @) and ratio rQ/(^{x) = 
bQ{x) /bc{x) from the hyperspherical treatment (see formulas (|34|) and (|35|) ) . as a function of the 
mass ratio x. 

FIG. 2: Ratio Rq/y{x) from the geometrical treatment (see formula and ratio rc/y(x) = 

bQ{x) /byix) from the hyperspherical treatment (see formulas H34|) and ()41|)). as a function of the 
mass ratio x. 

FIG. 3: Ratio Rm/y{x) from the geometrical treatment (see formula and ratio ?"M/y(^) from 
the hyperspherical treatment (see formula (|42|) '). as a function of the mass ratio x. 

FIG. 4: Integration domain of quantity L(R,6,x) ^ & function of angles 9 and x (see Sect. IIIl'X)l . 

FIG. 5: Ratio Ry/c ~ 1.086 from the geometrical treatment (see formula (jH))) and ratio ry/cix) = 
byix) /bc{x) from the hyperspherical treatment (see formulas (|4H) and (|35|) ). as a function of the 
mass ratio x. 
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